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Abstract 

Background: Pyrethroid resistance has been well documented in Anopheles arabiensis, one of the major African 
malaria vectors, and the predominant malaria vector in South Africa. 

Methods: In this study, the genetic basis of pyrethroid resistance in a selected laboratory strain of An. arabiensis 
from South Africa was investigated using a custom-made microarray, known as the An. gambiae detoxification chip. 

Results: A large number of P450 genes were over-transcribed, as well as a suite of redox genes and glutathione 
S-transferases. The five genes that showed the highest level of gene transcription when compared with an 
insecticide susceptible strain were: CYP6AG2, CYPZ1, TPX2, CYPZ2 and CYP6P1. 

Conclusions: Permethrin resistance in South African An. arabiensis is associated with increased transcription of 
multiple genes, and a large proportion of these genes were also previously recorded as over-transcribed in 
another An. arabiensis strain selected for resistance to DDT with cross-resistance to deltamethrin. The 
deltamethrin resistance developed de novo in the DDT-selected strain and is most likely due to increased 
transcription of those genes associated with DDT resistance. However, of particular interest was the fact that the 
strain selected for resistance to pyrethroids did not develop de novo resistance to DDT. These differences are 
compared and discussed. 

Keywords: South Africa, Insecticide resistance, Cytochrome P450, Glutathione S-transferase, Permethrin, 
Deltamethrin, Microarrays 



Background 

Pyrethroids are the most commonly used insecticides for 
malaria vector control They are most often used for in- 
door residual spraying, and are the only class of insecti- 
cide approved for treatment of bednets due to their low 
levels of mammalian toxicity [1-3]. However, pyrethroid 
resistance in African malaria vectors, including Anopheles 
arabiensis, is well documented and has been reported in 
many African countries. 

The Global Plan for Insecticide Resistance Manage- 
ment in Malaria Vectors (GPIRM) [1] provides a generic 
guideline for managing insecticide resistance where it 
occurs. Strategies for preventing and managing insecticide 
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resistance are imperative to ensure that the limited num- 
ber of insecticides available for vector control are pro- 
tected. As part of this strategy, insecticide resistance 
should be monitored through routine surveillance and the 
mechanisms of insecticide resistance should be charac- 
terised where possible [1]. Information on insecticide re- 
sistance mechanisms in malaria vectors can be obtained 
using several techniques that vary in sophistication and 
cost. The more costly methodologies are not feasible in 
many African countries where the disease burden is the 
greatest. 

Resistance to pyrethroids is mainly due to two mecha- 
nisms; enhanced enzyme detoxification or reduced target 
site sensitivity. Cytochrome oxidases (P450s), glutathione 
S-transferases (GSTs) and esterases are generally associ- 
ated with enzymatic detoxification of insecticides, while 
mutations in the voltage-gated sodium channel (commonly 
referred to as knockdown resistance or "kdf) effectively 
decrease target site sensitivity [4,5]. The occurrence of kdr 
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mutations is easily detected through the use of a polymer- 
ase chain reaction (PCR) assay [5] but this mechanism is 
not clearly correlated with all instances of pyrethroid re- 
sistance, especially deltamethrin resistance [6]. Identifying 
metabolic based resistance is more challenging and costly, 
especially when using microarrays [7]. The use of micro- 
array technology is particularly helpful as it allows for 
rapid analysis of numerous genes simultaneously, and as 
such provides substantial genetic information relating to a 
particular phenotype. 

Genes associated with pyrethroid resistance include the 
esterase, GST and P450 super families [8-13]. Metabolic 
detoxification of pyrethroids has been demonstrated in 
An. gambiae [12,14,15], An. arabiensis [8,10,16] and An. 
funestus (reviewed by [17]). Anopheles arabiensis is a 
major vector of malaria in Africa including South Africa 
[18,19]. Historically, this species was considered to be sus- 
ceptible to all classes of insecticide. In 1996 an extensive 
survey of insecticide susceptibility in An. arabiensis was 
conducted in three areas in South Africa (KwaZulu-Natal 
[KZN], Limpopo Province and the Kruger National Park) 
[20]. Bioassays revealed that samples of An. arabiensis 
from all three areas were susceptible to DDT and 
deltamethrin [20]. 

Current malaria vector control in South Africa is based 
on a mosaic approach using DDT and pyrethroids [19]. 
These interventions are the mainstay of mosquito vector 
control in South Africa even though An. arabiensis is both 
exophilic and endophilic. In 2003, An. arabiensis sampled 
from two areas in KZN were found to be resistant to 
DDT, but susceptible to deltamethrin [21]. 

Subsequent studies were initiated to investigate the 
molecular basis for the reported DDT resistance in KZN 
An. arabiensis. The use of laboratory strains to study re- 
sistance mechanisms carries the advantage of excluding 
factors that can confound the analysis (e.g. exposure to 
agricultural insecticides, effect of temperature, larval diet, 
etc). Laboratory based DDT selection pressure was exerted 
on this population over a period of several years and 
resulted in the fixation of the L1014F kdr mutation [8]. In 
addition to target site resistance, metabolic detoxification 
was also identified and numerous gene transcripts were 
over-expressed and associated with the resistant pheno- 
type [8]. 

Three years after the discovery of DDT resistance, pyr- 
ethroid resistance was also reported in wild caught An. 
arabiensis from KZN [16]. This population was colonised 
in 2005 and selected for resistance to permethrin in order 
to characterise the resistance mechanisms involved. Pyreth- 
roid resistance in this population (named KWAG-perm) 
was primarily associated with increased activity of mono- 
oxygenases (P450s) [16]. Analysis of six specific P450 
genes showed increased transcription for three genes in 
association with the pyrethroid resistant phenotype [22]. 



The aim of this study was to evaluate the levels of 
transcription of detoxification enzymes associated with 
the South African permethrin resistant strain (KWAG- 
perm) using a high throughput microarray analysis to 
supplement the information obtained by Munhenga and 
Koekemoer [22]. In addition to this, the transcript pro- 
files between the pyrethroid-selected (KWAG-perm, per- 
methrin resistant) and DDT-selected (MBN-DDT, DDT 
and pyrethroid resistance) strains were compared to iden- 
tify both strain specific transcripts as well as transcripts 
shared between these two strains. The implications of 
these results for malaria control are discussed. 

Methods 

Mosquito strains 

Pyrethroid resistant strain (KWAG-perm) 

The An. arabiensis permethrin resistant strain, KWAG- 
perm, and the equivalent susceptible strain, KW AG -base, 
were used for this experiment. KW AG -base originates 
from Mamfene in KwaZulu-Natal (KZN) province, South 
Africa. This population was characterised as permethrin 
resistant at the time of colonisation [16]. However, in the 
absence of selection by permethrin exposure, the labora- 
tory strain reverted back to being susceptible (after a 
period of one year). When permethrin selection was 
reintroduced, resistance rapidly returned (within four 
generations) and this strain was named KWAG-perm 
(Table 1) [16]. Mosquitoes were reared under standard in- 
sectary conditions (26°C ± 2°C, 70-80% relative humidity, 
and 12:12 light:dark cycle). 

DDT resistant strain (MBN-DDT) 

Details on this strain can be found in Nardini et al. [8]. 
Briefly, MBN-base was colonised from Mamfene, KZN, 
and a sub-strain was selected for resistance to DDT. The 
selected strain, MBN-DDT, is resistant to DDT (0% mor- 
tality) as well as permethrin (4%), deltamethrin (34%), 
bendiocarb (77.5%) and propoxur (65.3%) through cross- 
resistance mechanisms (Table 1). The strain is homozy- 
gous for the L1014F kdr mutation which is mainly respon- 
sible for the observed DDT and permethrin resistance [8]. 
The deltamethrin resistance is mainly due to enhanced 
monooxygenase enzymatic detoxification [8]. 



Table 1 Summary of the An. arabiensis strains used in the 
study 



Strain 


Selected on: 


Resistant to: 


KWAG-base 


No selection 




KWAG-perm 


Permethrin (0.75%) 


Permethrin 


MBN-base 


No selection 


DDT (low-level) 


MBN-DDT 


DDT (4%) 


DDT, deltamethrin and permethrin 
(cross-resistance) 
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RNA extractions and preparation of amplified mRNA 
Pyrethroid resistant strain (KWAG-perm) 

RNA extractions were carried out according to Nardini et al. 
[8]. Briefly, female mosquitoes were maintained on a 10% 
sugar solution and RNA was extracted from 15 mosqui- 
toes three days after emergence, representing one bio- 
logical repeat. RNA was extracted for three biological 
repeats using the PicoPure® RNA Isolation Kit (Arcturus) 
according to supplied methodology (with the inclusion of 
a DNase treatment, as described in the user manual). 
RNA was quantified using a NanoDrop and quality was 
assessed using an Agilent 2100 Bioanalyzer. Extracted 
RNA was stored at -70°C. In order to obtain larger quan- 
tities of mRNA for microarrays, the RiboAmp™ RNA 
Amplification Kit (Arcturus) was used. Amplified mRNA 
was quantified using a NanoDrop. 

Preparation of microarrays and data analyses 
Pyrethroid resistant strain (KWAG-perm) 

Three independent biological repeats were used, and for 
each biological repeat, two technical repeats were per- 
formed that included dye swaps in order to compensate 
for dye bias. Microarrays were prepared as described in 
Nardini et al [8]. 

The arrays were scanned using the Genepix 4000B 
scanner (Molecular Devices, USA). The photomultiplier 
tube (PMT) values were adjusted to give a pixel ratio of 
approximately 1. Spot quality and background intensities 
were examined and corrected using Genepix Pro 6.0 
software (Axon Instruments, USA). Saturated features 
were excluded from the analysis. 

Gene expression data were analysed using Limma 
version 2.12.0 (Bioconductor) [23] in R, version 2.8.0 
(http://cran.r-project.Org/bin/windows/base/old/2.8.0/). 
Raw intensity values for each spot were calculated, and 
then background corrected by the method "normexp" 
with an offset of 50 [24]. The corrected intensity values 
were transformed to log-ratios and then normalised. 
Composite Loess was used for within array normalisation. 
In this method, control spots and features, per sub-array, 
are used for producing non-linear, best-fit lines [25]. The 
use of control spots ensures that the resulting best-fit line 
is not biased by differential expression of genes. Con- 
versely, the use of all genes for normalisation improves 
stability with respect to the number of spots and, most 
importantly, provides flexibility in terms of print-tip group 
trends that might be observed where sub-array Loess 
curves are used [26] . Print- tip peculiarities were present in 
some slides, hence the choice of the normalisation 
method. The "Aquantile" method was used for between 
array normalisation. Genes with adjusted ^-values < 0.05 
and fold-changes (FCs) > 2.0 were considered to be statis- 
tically significant. These data have been deposited into 
Vectorbase (https://www.vectorbase.org). 



Quantitative real-time PCR (qPCR) 
Pyrethroid resistant strain (KWAG-perm) 

Real-time PCR was carried out in order to validate the 
outcome of the microarray experiments. RNA was 
extracted from 15 three day old An. arabiensis females 
(represents one biological repeat, and three biological re- 
peats were prepared) that had been supplied with 10% 
sugar solution. TRI-Reagenf solution ( Sigma- Aldrich) 
was used according to the supplied methodology with a 
DNase treatment included. Samples were reverse-transcribed 
into cDNA using the QuantiTecf Reverse Transcription 
Kit (Qiagen). The cDNA was stored at -70°C until required 
for qPCR. 

Two genes, CYP6AG2 and TPX2, were evaluated using 
real-time PCR. Beacon Designer™ (Premier Biosoft) or 
Invitrogen s free online primer design tool, OligoPerfect™ 
Designer, were used to design primers. These were based 
on An. gambiae sequence information (CYP6AG2: GB 
AY745225; TPX2: TIGR TC48596). The reference gene 
used for qPCR was based on that of Munhenga and 
Koekemoer [22], where studies were conducted in the 
laboratory on the same mosquito strains. These authors 
reported that 18S rRNA showed the most stable expres- 
sion of the six potential reference genes tested. The for- 
ward and reverse primer sequences for 18S were the 
same as those of Munhenga and Koekemoer [22]. Each 
PCR reaction comprised 12.5 \A IQ™ SYBR super-mix 
(Bio-Rad), 4 \A primer (primer concentration optimised 
for each primer set), 1 \A cDNA (100 ng/|il) and nuclease- 
free water to a final volume of 25 \A. Cycling conditions 
were as follows: 93°C for 3 minutes, followed by 35 cycles 
of (94°C/20 seconds, x°C/25 seconds, 72°C/30 seconds), 
and a final extension at 72°C for 10 minutes (followed by 
melt curve analysis). Primer sequences and annealing tem- 
peratures are described in Table 2. Data were analysed 
using the Pfaffl [27] method. For each gene of interest, the 
PCR product was sequenced in both directions, so that 
the presence of the correct product was confirmed (in 
addition to melt-curve analysis). Where necessary (for 
small amplicons i.e. < 110 base pairs), samples were 
cloned and then sequenced. 

Results 

KWAG-base and KWAG-perm have been evaluated for 
insecticide resistance on an ongoing basis in the labora- 
tory. The base strain showed very little resistance to 
permethrin (97.8% mean mortality 24 hours post expos- 
ure to 0.75% permethrin), while the resistant strain 
showed a high frequency of survivors on permethrin 
(42% mean mortality 24 hours post exposure to 0.75% 
permethrin) [22]. 

Microarray experiments indicated that 29 genes were 
over-transcribed according to the criteria outlined above. 
Most of these were P450 genes (55%), followed by redox 
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Table 2 KWAG-base/KWAG-perm primer information for 
qPCR (F = forward, R = reverse) 



Gene Primer sequence Primer Annealing 

concentration temperature 



CYP6AG2 


F 5'- TTG TGC TGC CGT 
ACT CG-3' 

R 5'- TAC TAT CGC CCG 
TCT CAC CT -3' 


2.0 pM 


59.4°C 


TPX2 


F5'- GGATGTTOTGG 
GGA ATA CG -3' 

R 5'- TGT GCG A^ AGC 
CTC CTC TT-3' 


3.5 mM 


56.3°C 


18S 


F5'-TAC CTGGGC G^ 
CTA CTC -3' 

R5'- OTTGA GCACTC 
TAA TTT G^ C -3' 


a 


b 



a Primer concentration was the same as that of CYP6AG2 and TPX4 when used 
as a reference gene in each case. 

b Annealing temperature was the same as that of CYP6AG2 and TPX4 when 
used as a reference gene for each. 



genes (21%), and GSTs (14%). A group of additional 
genes, such as cytochrome c, a ribosomal gene (RPS26) 
and a receptor protein (GPR npy 3) accounted for 10% 
of the over-transcribed genes. The four genes with the 
highest transcript abundance were CYP6AG2, CYP6Z1, 
TPX2 and CYP6Z2 (Table 3), in order of statistical sig- 
nificance (i.e. adjusted ^-value). 

Two significantly over- transcribed genes, CYP6AG2 and 
TPX2, were selected for qPCR validation of microarray 
data. The two other "top" genes (CYP6Z1 and CYP6Z2) 
were already evaluated with qPCR [22]. CYP6AG2 was the 
most significantly over-transcribed gene, and TPX2 is of 
interest because the TPX genes have been found to be im- 
portant in other studies reporting on insecticide resistance 
[8,10]. The genes CYP6AG2 and TPX2 produced fold 
change values of 2.9 and 4.5 respectively (Figure 1). The 
FC value of CYP6AG2 and TPX2 were comparable to 
those obtained by microarray analysis. 

Munhenga and Koekemoer [22] used qPCR to evaluate 
transcription of specific genes (identified through pub- 
lished literature) that have been implicated in permeth- 
rin resistance in An. arabiensis. These included CYP6Z1, 
CYP6Z2, CYP6Z3, CYP6M2, CYP6P3 and CYP4G16. Of 
these, five appeared on the list of over-transcribed genes 
in this microarray study (Table 3) and three of these 
genes [CYP6Z1 (4.7-fold), CYP6Z2 (1.7-fold) and CY- 
P6M2 (2.2-fold)] were significantly over-transcribed in 
the qPCR study [22]. The FC values of all six genes 
analysed by Munhenga and Koekemoer [22] were com- 
pared with the FC values obtained here for the same 
genes using microarray analysis (Figure 2). The FC 
values obtained by microarray analyses were compared 
with those obtained using qPCR for each gene using a 
£-test. The FC values of only three of the genes were 



significantly different to the FC values obtained by 
microarray analysis. These were CYP6Z2 (p = 0.0251), 
CYP6Z3 (p = 0.0219) and CYP6M2 (p = 0.0092). 

Transcript profile comparison between the DDT-selected 
strain and the pyrethroid-selected strain (this study) 

Both MBN-DDT and KWAG-perm originate from 
Mamfene (KZN). However, they were placed under dif- 
ferent insecticide selection pressures (depending on the 
insecticide resistance phenotypes present at the time of 
colonisation). The list of over- transcribed probes ob- 
tained in the permethrin-resistant strain (KWAG-perm, 
pyrethroid resistance only) was compared with that 
obtained in the DDT selected strain (MBN-DDT, pyreth- 
roid and DDT resistance [8]) (Figure 3). In MBN-DDT, 
kdr was closely associated with DDT and permethrin re- 
sistance, while deltamethrin resistance was largely based 
on enhanced enzymatic detoxification. Five GST transcripts 
were unique in the DDT-selected strain (MBN-DDT), while 
14 transcripts were unique to the permethrin-selected 
strain (KWAG-perm). Fifteen gene transcripts were 
shared between these of which 10 (66%) belonged to the 
P450-enzyme family. 

Discussion 

The use of controlled insecticide selection in the labora- 
tory aids in the study of changes in resistance status and 
its effect on certain adaptive features [28]. A list of 29 
over-transcribed genes according to the criteria used 
(adjusted ^-value < 0.05; FC > 2) was produced for the 
KWAG-perm strain. Most of the genes that were over- 
transcribed belong to the cytochrome P450 superfamily, 
supporting the findings of previous studies [16,22]. The 
four most statistically significant P450 genes included 
CYP6AG2, CYP6Z1, CYP6Z2 and CYP6P1, all of which 
belong to the CYP6 family and are well-known for their 
role in insecticide resistance [14,29]. The most signifi- 
cant gene, CYP6AG2, is not well characterised and it 
would be valuable to determine the functionality of this 
enzyme. This gene is located on the chromosome arm 
2R, a region that is not typically linked to pyrethroid re- 
sistance. Another gene in this family, CYP6AG1, was 
identified in a wild deltamethrin resistant An. arabiensis 
population from Cameroon [10]. CYP6Z1, CYP6Z2 and 
CYP6M2 showed increased transcription in this study as 
well as in the study by Munhenga and Koekemoer [22] . 

Pyrethroids are also known to cause oxidative stress by 
inducing lipid peroxidation, protein oxidation, and deple- 
tion of reduced glutathione [30] . This effectively increases 
the toxicity of the insecticide. A system of enzymes, in- 
cluding the superoxide dismutases (SODs), catalases, per- 
oxidases, cytochrome C and GSTs, is present and provides 
defense against these reactive oxygen species (ROS) 
[31,32]. If ROS are not metabolised, they damage 
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Table 3 List of probes over-transcribed in the resistant phenotype, KWAG-perm 



Gene 


Function 


Adj. p-value 


FC 


Accession number 


Location 


CYP6AG2 


Cytochrome P450 


2.58E-7 


4.1 


AY745224 


2R 


# CYP6Z1 (oligo) 


Cytochrome P450 


3.05E-7 


4.7 


AF487535 


3R 


TPX2 


Thioredoxin peroxidase 


4.56E-7 


2.3 


TIGR: TC48596 


3 L 


*CYP6Z2 


Cytochrome P450 


1.14E-6 


3.6 


XM_3 17252 


3R 


CYP6P1 


Cytochrome P450 


1.14E-6 


2.2 


AY028785 


2R 


GSTU1 


Glutathione S-transferase 


3.03E-6 


2.2 


XM_309135 


X 


S0D2 


Superoxide dismutase 


3.81 E-6 


3.5 


AY524130 


2 L 


CYP12F2 


Cytochrome P450 


3.81E-6 


2.9 


AY1 76050 


3R 


CYP6Y2 


Cytochrome P450 


3.81 E-6 


2.2 


AY1 93728 


3R 


GPR npy 3 


G protein coupled receptor 


4.83E-6 


3.2 


ENSANG: 


2R 










G00000009317 




CYP9J5 


Cytochrome P450 


9.14E-6 


3.2 


AY748830 


3 L 


GPX1 


Glutathione peroxidase 


9.14E-6 


2.4 


AY842257 


2R 


*CYP6P3 


Cytochrome P450 


1 .68E-5 


2.6 


AF487534 


2R 


# CYP6Z1 (cDNA) 


Cytochrome P450 


1 .96E-5 


3.0 


AF487535 


3R 


CYP6M3 


Cytochrome P450 


1 .99E-5 


4.5 


AY1 93730 


3R 


S0D1 


Superoxide dismutase 


2.35E-5 


2.7 


AY505417 


3 L 


*CYP6Z3 


Cytochrome P450 


2.95 E-5 


3.0 


AY1 93727 


3R 


CYP6AK1 


Cytochrome P450 


5.42E-5 


3.6 


AY745227 


3 L 


CYP9M1 


Cytochrome P450 


1.03E-4 


2.6 


AY748836 


3R 


CYP12F4 


Cytochrome P450 


1 .68E-4 


2.3 


AY1 76048 


3R 


CYP4H24 


Cytochrome P450 


2.05 E-4 


2.6 


AY062206 


X 


*CYP6M2 


Cytochrome P450 


3.38E-4 


4.6 


AY1 93729 


3R 


GSTD1-3 


Glutathione S-transferase 


6.83 E-4 


2.0 


AF071163 


2R 


CYP6M 1 


Cytochrome P450 


1.16E-3 


2.2 


AY062208 


3R 


GSTS1-1 


Glutathione S-transferase 


1 .64E-3 


2.5 


L07880 


3 L 


RPS26 


Ribosomal protein 


2.16E-3 


3.1 


EMBL: 


3 L 










4A3A-AAM-G-11-R 




GPX3 


Glutathione peroxidase 


2.16E-3 


2.1 


AY745228 


X 


TPX4 


Thioredoxin-dependent peroxidase 


3.60E-3 


2.8 


AY745235 


3 L 


GSTS1-2 


Glutathione S-transferase 


4.15E-3 


2.3 


AF5 13639 


3 L 


Cytochrome_C 


Cytochrome c 


2.17E-2 


2.0 


TIGR: TC48590 


3R 



Probes are listed in order of significance (adjusted [adj.] p-value), and then by fold change value. The accession number refers to Genbank, unless otherwise 

specified in the table; FC Fold change, E Exponent. 

# Gene transcript validated by qPCR by Munhenga and Koekemoer [22]. 



important compounds such as lipids, proteins, nucleic 
acids and carbohydrates, and ultimately cause cell death 
[32]. A number of redox enzymes showed elevated levels 
of transcription in KWAG-perm, which may be associated 
with resistance to pyrethroids. These included thioredoxin 
peroxidases (TPX2 and TPX4), glutathione peroxidases 
(GPX1 and GPX3), superoxide dismutases (SOD1 and 
SOD2) and cytochrome C Superoxide radicals are con- 
verted to hydrogen peroxide and oxygen by SODs, and 
hydrogen peroxide is converted to water and oxygen by 
catalases; or to water, by peroxidases [33] . 



Four GSTs {GSTU1, GSTD1-3, GSTS1-1, GSTS1-2) 
were also over- transcribed in KWAG-perm. Over- 
transcription of GSTs has been observed before in pyreth- 
roid resistant An. arabiensis [10], An. gambiae [7,9,34] 
and in other insects [35-37] . In conjunction with glutathi- 
one, GSTs function as antioxidants by limiting peroxida- 
tion and by limiting (termination) "free-radical cascades" 
(see [30]). In addition, the GSTs are able to bind to pyre- 
throids and provide protection by sequestration [37] . 

In the recent study of Nardini et al. [8], the authors 
identified genes that were over-transcribed in a DDT 
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Figure 1 The volcano plot of KWAG-perm versus KWAG-base microarray data. The cut-offs for significance are indicated by the dotted lines 
(adj. p-value < 0.05; FC > 2) and the top 6 genes that met this criteria have been labeled. All positive FC values belong to the genes that are 
over-transcribed in the resistant strain (KWAG-perm), while all negative FC values belong to genes that are over-transcribed in the 
susceptible strain. 



and pyrethroid resistant strain (MBN-DDT). These 
genes were mainly found to be associated with the ob- 
served cross resistance to deltamethrin (the use of syner- 
gists clarified the roles of the over- transcribed genes, 
and implicated kdr in DDT and permethrin resistance). 
The gene transcripts that were identified and which are 
shared between the permethrin-selected (KWAG-perm) 
and DDT-selected strains (MBN-DDT), might explain 
the protection against deltamethrin in the DDT selected 
strain. In the process of selecting for DDT resistance it 
appears that there was also selection for increased tran- 
scription of genes (n = 15) that coincidentally protect 
against pyrethroids as well. However, the converse was 
not true and the pyrethroid selected strain did not de- 
velop DDT resistance. This implies that DDT selection 
pressure can also result in the de novo development of 
resistance to pyrethroids. Very few reports of DDT re- 
sistance without resistance to pyrethroids are available 
(e.g. [38]), which supports this hypothesis. 

It is possible that the genes that are unique in each 
instance (KWAG-perm versus MBN-DDT) are more 



important for resistance to either type I (e.g. GSTD2) or 
type II pyrethroid (CYPZ1), while those that are com- 
mon are likely to play a role in resistance to both kinds 
of pyrethroid (e.g. CYP6M2). Those genes that are 
shared between the pyrethroids and DDT selected strain 
might be able to metabolise both classes of insecticides 
(DDT as well as pyrethroids). This is specifically true for 
CYP6M2 [39], but to our knowledge similar metabolic 
studies are not available for the other 10 genes. The role 
of genes conferring carbamate resistance identified in the 
DDT-selected strain still needs to be clarified [8]. How- 
ever, P450s have been shown to provide protection against 
both pyrethroids and carbamates in An. funestus [40] . 

Conclusions 

The challenge with detoxification based resistant mecha- 
nisms is that it is not possible to develop a single mo- 
lecular assay for the detection of resistance (as in the 
case of target-site mediated resistance, for example, kdr 
detection). Furthermore, in target-site resistance mecha- 
nisms, genotype usually correlates with phenotypic out- 
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Figure 2 A comparison between average fold change (FC) +SEM values recorded using microarrays and qPCR. qPCR FC values included 
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values were compared to each other by means of a f-test. Pairs of bars marked with a * are significantly different. 
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over-transcribed in the deltamethrin-resistant strain. Five of the over-transcribed genes were unique to the deltamethrin strain and all of these 
were GSTs. 
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come. In the case of enzymatic resistance, this is not neces- 
sarily true. For example, three GSTs were over- transcribed 
in the strain here (the activity of GSTs is known to confer 
DDT resistance in mosquitoes [41,42]), but no DDT re- 
sistance is present in KWAG-perm. In addition, different 
genes were over-transcribed in the permethrin and delta- 
methrin resistant strains. 

Microarrays are extremely useful for identifying a large 
number of genes that are associated with a particular re- 
sistance phenotype and in particular, transcripts that are 
unique to a particular phenotype. However, this techno- 
logy is extremely expensive and not available in most 
African countries and an alternative to this approach is 
needed. Due to this problem, the results from the costly 
microarray approach used in this study were compared 
against those obtained by a less expensive qPCR study 
using the same An. arabiensis strain [22]. Selection of 
genes to be included in this qPCR study were identified 
through published literature [7,9,10,14,34]. This approach 
was used successfully by Munhenga and Koekemoer [22] 
who found >1.5 FC in transcription in four of the six 
genes that were selected for analysis. Five of the six genes 
from Munhenga and Koekemoer [22] were also identified 
through this microarray study. However, CYP6AG2 and 
TPX2 were novel in this study. 

These data suggest that relatively inexpensive qPCR 
studies can successfully be used to identify increased tran- 
scription of specific metabolic genes, especially in coun- 
tries that do not have access to advanced molecular 
systems eg. microarray analysis. However, many transcripts 
might be overlooked using qPCR due to the specificity of 
the technique. Microarray studies will in future expand the 
number of transcripts available for qPCR screening. 

The GPIRM has initiated a "five pillar" approach against 
insecticide resistance. One of these pillars is to gain know- 
ledge on mechanisms of insecticide resistance and on the 
effects of current insecticide resistance management strat- 
egies [1]. The present study and similar studies, contribute 
to achieving this goal. 
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